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It is shown that Schrodinger equation with combination of three potentials U — ~ar~^ + Pr + kr^ , 
Coulomb, linear and harmonic, the potential often used to describe quarkonium, is reduced to 
I a bi-confluent Heun differential equation. The method to construct its solutions in the form of 

^SJ . polynomials is developed, however with additional constraints in four parameters of the model, 

a, (3, k, I. The energy spectrum looks as a modified combination of oscillator and Coulomb parts. 
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I. INTRODUCTION 

Gauge theories of the strong interactions suggest that the coupling between quarks is weak at short distances 
r*! , but becomes very strong at large distances. This 'explains' the paradox where quarks appear to behave as quasi-free 
particles within hadrons but cannot be liberated from the hadrons. On the basis of these arguments it was conjectured 
that heavy quarks would move nonrelativistically within hadrons. Thus bound states of a heavy quark and antiquark 
should be a hadronic analogue of the positronium system of a bound electron and positron. 
I This so-called quarkonium system might then be interpreted according to the familiar rules of nonrelativistic 
^-H ■ quantum mechanics using a potential to describe the interquark force (see and references therein). 
J> ' At very small distances, the potential is expected to take a form like the Coulomb force, corresponding to the 
, exchange of a single massless gluon. At very large distances, a linear confining potential or something more complex 
seem to be appropriate. We will use a combination of three well-known potentials 

U = -- + /3r + kr'^ 

O ' ^ 

I ' that has some advantages for analytical treatment of the quarkonium problem in terms of Heun functions (this class 
y—i . of special functions b eing next extension to hypergeometric functions has become of primary importance in many 
' physical problems [3j-[32|). 

H II. SCHRODINGER EQUATION 

In Minkowski space, parameterized by spherical coordinate 

dS^ = c^dt^ - dr^ - r^{d9^ + sin^ Odcf)'^) , 
the Schrodinger equation for an arbitrary spherical potential has the form 



■^9 ^ 



(1) 



After separation of the variables, ^ — e 'i Yim{0, (f) R{r) , we get 

h? (I d^^d + 



2M \r^dr^ dr ' ^ ^^""^ 



R. (2) 



For a potential specified by (we assume the Coulomb attraction, so a > 0) 

a 



U = Pr + kr^ 

r 



eq. ^ reads 



d2 2d III + 1) ^ 2M 



2M f a ^ , 

+ ^ -B H I3r-kr- 



R = 



We may simplify notation by using (instead of more short form 



(f 2d 



— - + 2e+-- - pr-kr'^\R^Q 



and further with the substitution R{r) = r ^ f{r) we get 



d^ a 1(1 + 1) ^ , 2 , , 

— + 2e+ -Pr-kr^] f = 

ar'' r 



Let us see behavior of the curve 



P^(r) = 2e+^-^^-/3r-fcr^ 



To proceed further, let us examine the classical turning points - the roots of the equation 

- kr^ - I3r^ + 2tr^ + ar - /(/ + 1) = ~k{r - ri){r ~ T2){r - r3)(r - r4) = 
From ([71) it follows 

^ ^ ^ ^ 

- r = ''1 + ^2 + f 3 + ?'4 , 

k 

2e 

--^ = rir2 + J'lTa + rir^ + r2rz + r2r4 + r3r4 , 
a 

- = r2r3r4 + rirj,ri + rir2ri + rirara , 



k 



There exist the possibility when two roots are negative and two other are positive (see Fig. ] 
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FIG. 1: Finite classical motion at r G [?'3,?'4] 



Prom (O it follows two systems: 



^ = rir2 + ri(r3 + + r2(r3 + + r3r4 , 
k 

- = n r3r4 + r2 r3r4 + rir2(r4 + r3) , 



and 



r\r2 



f 3 - '"4 
fc r3r4 



Making the change of variables y — iKr (K ~ k^l'^): 



and using the substitution 



d^R 2dR (2e la Ijl + 1) i(3y 2,o^n 
dy^ y dy Ky ^ ' 



R = y^e^ye^y'F{y), 



we get 



d'^F 
dy"^ 



2(A+1) 

y 

2e 



ACy + 2B 



dF 
dy 



{AC^-l)y^+[4BC--^ 



+AAC + - — +6C 
K ^ 

With the choice of special values 



+ A-l{l + l) 2ABK + 2BK ~ia 



y 



yK 



F 



A = +l,-{l + l), C^±\, B^±^, 



take those related to possible bound states 



A = +1, / - r' 



C 



B 

eq. (jlip becomes simpler 

+ AAC + B'^ 
It is convenient to turn to the variable 



1 

= +2' 
2K^ 



<PF_ 
dy^ 



2(A+1) 



y 



ACy + 2B 



2t ^ 2ABK + 2BK - ia 

+ + Vk 



y 



dF 
dy 

F^O. 



z = — = Kr 
i 



then eq. (|T4l) reads (remember that A = +1 C = +1/2 

d^F 
dz^ 



+ [-2-{2l + l) + {^-B^) 



-2z + 2.B+i±i^)^ 
z J dz 

-{-2iB){l + l) + a/K 



F = 0, 



which can be recognized as a biconfluent Heun equation for functions H{a, h, c, d, z 
(fH ( l + a\ dH ( -h{a + l)/2- d/2 



dz^ \ z ) dz 

with parameters 



2z-6+^— —+ -2- a + c+ — ^ " —\H = Q (16) 



Let us present solutions of eq. p6|) as a series 

= 1 + CiZ + 02^2 + C3Z^ + •■• = X! ' 

n=0 

ff'(z) = Ci + 2C2Z^ + Scgz^ + ... = ^ nc„0""^ , 

Ji=l 

CXD 

77"(z) = 2x1 C2Z° + 3x2 cgz^ + ... = ^ n{n ~ l)c„z""2 . 



n=2 



Eq. (Uni) gives (let D = -b{a + l)/2 - d/2) 



d^H ( , l + a\ dH f D\ 

-2z-b+ — + { -2- a + c+ —\ H 



dz"^ \ z j dz 



z 



that is 



n{n - l)c„z" 



n=2 

+ (-2z - b + — ^) y nc„z"-i + 

oo 

+ (-2-a + c+-)y c„z" = 0. (18) 

n=0 



Prom (UHl) it follows 



y n(n - l)c„z""^ - y 2nc„z" - y 5nc„z""^ + 

Ti— 2 71—1 n—l 

oo oo oo 

+ y (1 + a)nc„z"-2 + ^(_2 - a + c)c„z" + y Dcr.z'^-^ 

n—l n— n— 

oc oo oo 

y (iV + 2){N + l)cAr+2^" - y 2NcNZ^ - y 6(A^ + l)cjv+iz^ + 

oo oc oo 

+ y (l + a)(iV + 2)cAr+2Z^+ y (-2-a + c)cArz^+ y Dcn+iz^ . 



N=-l N=0 N=-l 
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Collecting similar terms 

oo 

2c2 +^{N + 2){N + 1)CN+2Z^ - 

oo oo 
- J2 ^NcnZ^ - feci - ^ b{N + l)cAr+iZ^ + 
N=l JV=1 

oo 

+(1 + a)ciz-'^ + (1 + a)2c2 + ^ (1 + a){N + 2)cn+2Z^ + 

JV=1 

oo 

+(-2 - a + c)co + ^ (-2 - a + c)cnz^ + 

JV=1 

oo 

+Dcoz~'^ + £>ci + ^ Dcn+\z^ = 



JV=1 



we get 



[ (1 + a)ci + Dco] + 
[ 2c2 - feci + (1 + a)2c2 + (-2 - a + c)co + i'ci ] + 



+ ^ [ (A^ + 2)(7V + l)cjv+2 - 2NcN - b{N + l)cjv+i + 

JV=1 

+(1 + a){N + 2)cjv+2 + (-2 - a + c)civ + i'cjv+i ] -z^ = . 
Thus we arrive at the recurrent relations for series coefficients 

D (2 + a - c)co + (fe - £')ci 
(l + a) (o + 2)2 
(2jV + 2 + g - c) + [-i? + fe(jV + 1)] c^+i 
"^+^ = iN + 2){a + N + 2) ' ^=1'2'- (19) 

Prom this, after simple change in notation, we obtain 

D (2 + a - c)co + (fe - -D)ci 

(1 + a) [a + 2)2 

_ {2n + a- c) Cn-i + {-D + fen) c„ „ _ o q 4 ^90^ 

Cn+i 7 — — TTT — : — TT , n- 2,6,4:,... (2Uj 

(n + l)(a + n + 1) 



Remember that 



we obtain 



, , /3 , 2a 2e fe2 

a = 2l + l, b=j^, d=-— c=^ + - 



-D^h(l + l)-^, h-D = h(l + 2)-^, 
2b-D = b{l + 3)-^ , ... fen - D = fe(/ + n + 1) - ; (21) 

note that the energy parameter does not enter relations in (2.17), instead it is presented only in the parameter c 
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Recursive relation can be rewritten in the form 

D 



C2 



C3 



C4 



C5 







2(2/- 


r Cn 

h2) ° 


[(2a 


+ 1) + (1 - c) ] Co + (6 - 






2 [(2(^ + 1) - 






m + 


2) + (1 - c) ] Cl - 


+ (25- 






3 [2(? + l) + 


2] 




[2(3 4 


- + (1 - C) ] C2 - 


-f (35- 


D) C3 




4 [2(; + l)4 


■3] 




[2(4H 


- + (1 - C) ] C3 


+ (45- 


-D)C4 




5 [2(^ + 1) ^ 


-4] 





_ [ 2(n + + (1 - c) ] Cn-i + (n5 - D) c„ 

(n + 1) [2(/ + l) + n] - ^^^^ 

Ont can try to obtain polynomial solutions. The first simplest possibility can be realized by adding two constraints: 

ci = 0, C2 = (24) 

which gives 

f5{l + 1) 

c=(2(/ + l) + l =^ |i + ^ = 2(/ + l) + l. (25) 

Result (l25l) may be rewritten as follows 
the case rt = 

Let us consider the case n = 1 which is realized by adding two constraints 

C2 = , [ (2(/ + 1) + (1 - c) ] Co + (5 - D)c^ = , 

C3 = , 2(/ + 2) + (1 - c) = , (27) 



which leads to 



In explicit form they are 



and 



c = 2(/ + 2) + l, 
2+^^-^)2(2^-«- 



2e 52 

_ + _^2(/ + 2) + l, (29) 



4(2/ + 2) - [--| + 5(/ + 2)][--| + 5(/ + 1)] = . (30) 



Eq. (|3Q)) is a 2-order equation with respect to 6 

2 

52(/ + !)(/ + 2) - 5-^(2/ + 3) - 4(2/ + 2) + - 
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with its solution 



a a + 3/2) 
^K{l + l){l + 2) 



(; + 3/2) 
^/T^ V(^ + l)(^ + 2) 



+ 



4(2/ + 2) - aV.ft:2 
(Z + l)(Z + 2) 



1/2 



or 



1/4 



+ ■ 



(/ + !)(/ + 2) y^JC' (Z + l)2(/ + 2)2 {1 + 2)' 
Let us consider several first coefficients of the bi-confluent Heun series: 



D 



ci = - 



C2 



(2(/ + l) + (l-c) 
2 [ (2(/ + l) + l 



Co 



2(2/ + 2) 

{h-D) 



Co ; 



{2{l + 1) + (1 - c) 
2 [ (2(Z + 1) + 1 ] 



Co + 



{h-D) 



2 [(2(Z + 1) + 1 
D 



2[(2(Z + 1) + 1] V 2(2Z + 2) 



ci 



Co 



C3 



{2b - D) 



2(^ + 2) + (l-c) 
3 [2(; + l) + 2] ' 3 [2(/ + l) + 2] 
2(; + 2) + (l-c) / D 



C2 = 



+ 



{2b - D) 



3 [ 2(Z + 1) + 2 ] 



3 [2(; + l) + 2] 
(2(; + l) + (l-c) {b-D) 
2[(2(« + l) + l] °^ 2 [(2(Z + 1) + 1] 
2(3 + Z) + (l-c) 



Co 



Co 



C4 



2(3 + /) + (1 -c) 
4 [2(/ + l) + 3] 



{2{l + 1) + (1 - c) 



4 [2(/ + l) + 3] 

(7; - D ) 



C2 



2 [(2(/ + l) + l] 

(35 - D) 



Co 



2 [(2(/ + 1)4-1] 
2(l + 2) + (l-c) 



2(21 + 2) 
D 

2(2/ + 2) 
(3b - D) 
4 [2(/ + l) + 3] 
D 

Co 



C3 



4[2(/ + l) + 3] (3[2(Z + l) + 2] 



+ 



{2b - D) 



3 [2(Z + l) + 2] 



{2{l + 1) + (1 - c) 
2 [ {2{l + 1) + 1 ] 



Co + 



{b-D) 



2[(2(Z + 1) + 1] V 2(2Z + 2) 



2(2/ + 2) 
D 

"2(21 + 2) 
D 



Co 



Co 



(31) 



The coefficient c„ represents a n-polynomial with respect to parameter h. 

In principle, it is easily to extend the above polynomial-based approach to general case. Indeed, let it be 



which gives 



From whence it follows 



c„+i = 0, 



Cn+l = 



Cn+2 = , 



[ 2(n + + (1 - c) ] c„_i + {nb - D) = , 
c„+2 = 0, [2(n+l + /) + (l-c)]c„ = 0. 



(l-c) = -2(n+l + 0, 
-2c„_i + {nb-D)cn = 0. 



(32) 



(33) 



(34) 



The problem is reduced to rather complicated polynomials. This method provides us with the formula for energy 
levels; however we obtain some additional constraints for four parameters, a,P,k,l (in the form of n-polynomial). 
This means that we are able to construct solutions in the polynomial form, but only at some special values of a, 13, k, I. 
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General structure of the 3-terni recurrent relations can be presented in a more short notation 

Ci = AqCo , 
C2 = EqCq +A1C1, 

C3 = Eici + A2C2 , 
C4 = E2C2 + A3C3 , 
C5 = £'303 + A4C4 , 



where 



Cn+l — En-lCn-1 + AnCn 



(35) 



_ 2(n + + (l-c) 

(n + l)[2(/ + l) + n]' (n + 1) [2a + l) + n] ■ ^^""^ 
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